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Abstract 

We analyze the self-similarity approach applied to study the hadron production in pp 
and AA collisions. This approach allows us to describe rather well the ratio of the proton 
to anti-proton yields in A-A collisions as a function of the energy at a wide range from 
a few GeV to a few TeV. We suggest a modihcation of this approach to describe rather 
well the inclusive spectra of hadrons produced in pp collisions at different initial energies 
from the AGS to LHG. 


1 Introduction 

The description of hadron production using statistical models has been pioneered several 
decades ago by E.Fermi pQ, I.Pomeranchuk [2], L.D.Landau [3] and R.Hagedorn |1]. The 
transverse momentum spectrum of particles produced in hadron-hadron collisions can be 
presented in a simple form ph ~ exp(—m/it/T), where ruht is the transverse mass of the 
hadron h and T is sometimes called the thermal freeze-out temperature. 

As it is well known, the statistical (thermal) models have been applied successfully to 
describe hadronic yields produced in heavy-ion collisions (see, for example, 0 - 0 and 
references therein). The temperature obtained in these analyzes is often referred to as 
chemical freeze-out tempeTa.tnTe and is consistently slightly higher than the thermal freeze- 
out temperature. At the same time, the source of very fast thermalization is currently 
unknown and alternative or complementary possibilities to explain the thermal spectra 
are of much interest. 

There exists a rich and wide variety of distributions covering a large range of applica¬ 
tions [T0HI2]. Those having a power law behaviour have attracted considerable attention 
in physics in recent years but there is a a long history in other helds such as biology and 
economics |13] . 

In high energy physics the power law distributions have been applied in [T4IIT3] to the 
description of transverse momenta of secondary particles produced in pp collisions. Indeed 


1 


the available range of transverse momenta has expanded considerably with the advent of 
the Large Hadron Collider (LHC). Collider energies of 7-8 TeV are now available in pp 
collisions and transverse momenta of hundreds of GeV are now common. Applications of 
the Tsallis distribution to high energy e^e“ annihilation have been considered previously 
in [19] . A recent review of power laws in elementary and heavy-ion collisions can be found 
in ESI- The modification of the Tsallis distribution [^12^ and its successful application 
to the analysis of the LHC data on the multiple hadron production in pp collisions at the 
central rapidity region has been done recently in |25j . 

There are other approaches, like the quark gluon string model (QGSM) |2SII2Z] or 
the Monte Garlo (MG) versions of string model [281129] . which are applied to analyze the 
hadron production mainly at the non central rapidity region. 

The inclusive spectra of hadrons produced in central pp collisions at low and large 
hadron transverse momenta pt were analyzed within the modified quark gluon string 
model (QGSM) and the perturbative QGD (PQGD) [SU]. There was suggested a con¬ 
tribution of the nonperturbative gluons at low transfer momentum squared [SniEI], 
which results in the satisfactory description of LHG data on the hadron pt-spectra at the 
mid-rapidity region and different energies (\/s). The description of the energy depen¬ 
dence of these spectra is an advantage in comparison with many theoretical approaches 
mentioned above. 

Almost all theoretical approaches operate the relativistic invariant Mandelstam vari¬ 
ables s, f, u to analyze the hadron inclusive spectra in the mid-rapidity region. As usual, 
the spectra are presented in the factorized forms of two functions dependent of t or p'f 
and s. However, there is another approach to analyze multiple hadron production in pp 
and AA collisions at high energies, which operates the four velocities of the initial and 
hnal particles [32]. It is the so called the self-similarity approach, which demonstrates a 
similarity of inclusive spectra of hadrons produced in pp and AA collisions, as a function 
of similarity parameter H. In fact, this approach is valid not in the complete kinematical 
region. That will be discussed in our paper. The hadron inclusive spectra obtained within 
this approach are presented as a function of the relativistic invariant similarity parame¬ 
ter, which can be related to variables t and s. The general form of such spectrum is not 
factorized over t and s. That is an advantage of the similarity approach in comparison to 
all the models mentioned above. However, the unfactorized form of inclusive spectra is 
significant at not large initial energies and it becomes independent of ^/s at large ^/s like 
the ISR, SPS and LHG energies. 

The approach of studying relativistic nuclear interactions in the four velocity space 
proved to be very fruitful [33] • In this article, we present a further development of this 
approach. 

2 The parameter of self-similarity. 

Within the self-similarity approach [321133] the predictions on the ratios of particles pro¬ 
duced in AA collisions at high energies were given in [33] • Let us briefly present here the 
main idea of this study. Gonsider, for example, the production of hadrons 1, 2, etc. in the 
collision of a nucleus I with a nucleus H: 

1 + 11^1 + 2 + ... ( 1 ) 

According to this assumption more than one nucleon in the nucleus I can participate in 
the interaction ([T]). The value of Nj is the effective number of nucleons inside the nucleus 
/, participating in the interaction which is called the cumulative number. Its values lie 
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in the region of 0 < iV/ < A/ [Aj - atomic number of nucleus I). The cumulative area 
complies with Nj > 1. Of course, the same situation will be for the nucleus //, and one 
can enter the cumulative number of Nu. 

For reaction © with the production of the inclusive particle 1 we can write the 
conservation law of four-momentum in the following form: 

{NiPj + NnPn-pif = 

{NirriQ + Nnmo + Mf, ( 2 ) 

where Nj and Njj the number of nucleons involved in the interaction; Pj, Pjj , pi are 
four momenta of the nuclei I and II and particle 1, respectively; mo is the mass of the 
nucleon; M is the mass of the particle providing the conservation of the baryon number, 
strangeness, and other quantum numbers. 

In [35] the parameter of self-similarity is introduced, which allows one to describe the 
differential cross section of the yield of a large class of particles in relativistic nuclear 
collisions: 


n = mm[^[{uiNi + (3) 

where uj and un are four velocities of the nuclei / and II. The values Nj and Nu 
will be measurable, if we accept the hypothesis of minimum mass mo(MiA^i +U 2 N 2 )‘^ and 
consider the conservation law of 4-momentum. Thus, the procedure to determine Nj 
and Nu, and hence If, is the determination of the minimum of If on the basis of the 
conservation laws of energy-momentum. 

Then, it was suggested [3ll|35] that the inclusive spectrum of the produced particle 1 
in AA collision can be presented as the universal function dependent of the self-similarity 
parameter If, which was chosen, for example, as the exponential function: 

Ed^a/dp^ = 

■ exp(-n/C' 2 ), (4) 

where a{Ni) = 1/3 -|- Nj/S, 

Oi{Nu) = 1/3-1- Njj/ 3, 

Cl = 1.9 ■ 10‘^mb ■ GeV~‘^ ■ C • st~^ and 
Cs = 0.125 ±0.002. 

3 Self-similarity parameter in the central rapidity re¬ 
gion 

In the mid-rapidity region (y=0, y is the rapidity of particle 1) the analytical form for If 
was found in [3l|. In this case Nj and Njj are equal to each other: Nj = Njj = N. 

A^ = [1±(1±<F5/$2)^/2]<1), (5) 

where 

$ = 2mo{mtchY + M)/sh‘^Y, (6) 

sh'^Y). (7) 

Here rriit is the transverse mass of the particle 1, rriit = {'ml C- rapidity of 

interacting nuclei. 
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And then - 


U = N- chY. 


( 8 ) 


This formula is obtained by searching for the minimum of the expression for If by dif¬ 
ferentiation of the right-hand side of equation (3) on variables Nj and Nu taking into 
account the conservation law (2). 

For baryons we have 

Iff, = {uiitchY — mi)chY/ {mosh^Y) (9) 

and for antibaryons - 

ffa = {mitchY + mi)chY/{mQsh‘^Y). (10) 


The results of calculations for the ratio of the antiproton cross section to the proton one 
after integration of Eqs. flOllTOl) over dmu are presented in Fig. [TJ This ratio is compared 
to the experimental data including the latest data at LHC [361138] . 



Figure 1: The dependence of ratio of the antiproton cross section to the proton one as a 
function of initial rapidity Y or energy {\fS, GeV) of the interacting nuclei.The points 
are the experimental data [5B1155] . 


One can see a very good agreement of the results of our calculations with the experi¬ 
mental data in the wide energy range. 

The same calculations can be made for other antiparticles and particles. However, 
there are poor experimental data in the central rapidity region. 

4 Further development of Baldin’s approach. 

As it is mentioned above, the exponential form for the hadron inclusive spectrum given 
by Eq.(jl]) was chosen as an example and using it we can satisfactorily describe the ratio 
of total yields of antiprotons to protons produced in heavy nucleus-nucleus collisions. 
Unfortunately, this simple form Eq.(|l]) contradicts to the LHC data on the inclusive 
spectra of hadrons produced in the central pp collisions, as shown in [30]. Therefore, we 
use the results of [30] to present the inclusive relativistic invariant hadron spectrum at the 
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mid-rapidity region and at not large hadron transverse momenta pt in a more complicated 
form, which consists of two parts. The first one is due to the contribution of quarks, which 
was obtained within the QGSM [2^127] using the AGK (Abramovsky, Gribov, Kanchelli) 
cancellation [32] of n-pomeron exchanges for inclusive hadron spectra at the mid-rapidity 
region. It is written in the following form [30]: 

OO 

E{d^a/d^p)q = (i)q{y = 0,pt) • ^[na„(s)] = 

n=l 

0g(2/ = 0,Pt)^(s/So)^ , (11) 

where o'n(s) is the cross-section for production of the n-pomeron chain (or 2n quark- 
antiquark strings); g = 21 mb - constant, which is calculated within the ”quasi-eikonal” 
approximation mi; So = 1 GeV^] A = [Q;p(0) — 1] ~ 0.12, where Q;p(0) is the sub critical 
Pomeron intercept [2^1271110] . 

The second part of the hadron inclusive spectrum at the mid-rapidity region was 
introduced in [3D1EI] assuming the contribution of the nonperturbative gluons and calcu¬ 
lating it as the one-pomeron exchange between two nonperturbative gluons in the collided 
protons [31]. This part was written in the following form [30] : 

OO 

E{d^a/d^p)g = 4)g{y = 0,pt) ■ ^(n - l)cr„(s) = 

n=2 


( OO OO \ 

^nan(s) - ^a„(s) j = (j)g{y = 0,pt) ■ [g{s/ sq)^ - and] , (12) 

71=1 71=1 / 

where and is the non diffractive pp cross section. 

Thus, taking into account the quark and gluon contributions we will get the following 
form for the inclusive hadron spectrum: 

E{d^a/d^p) = [(j)q{y = 0,pt) + (j)g{y = 0,pt) ■ {1 - and/9{{s/so)‘^)] ■ g{s/so)^ (13) 


The question arises, what is a relation of the similarity parameter 11 to the relativistic 
invariant variables s,p^ ? This relation can be found from Eqs.([5]lH]) using ch{Y) = 
y^/(2mo). Then, we have the following form for 11: 


n 


\ 2mo5 ^/s5 j 



(14) 


where 5 = 1 — drup/s; mu = \/pt + is the transverse mass of the produced hadron h. 
At large initial energies y/s » 1 GeV the similarity parameter 11 becomes 


n 


mu 

2mo(l — 4mo/s) 



+ 


M2 


m 


m 


1(1 — 4mQ/s) 


It 


(15) 


For TT-mesons mi = is the pion mass and M = 0; for iP“-mesons mi = mx is the kaon 
mass and M = mx] for iP^-mesons mi = mx and M = m\ — mx, m\ is the mass of the 
A-baryon. For vr-mesons sX pi » mf we have: 


n ~ 


mu 

mo(l — Xml/s) 


(16) 
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One can see that in a general case the similarity parameter 11 depends on p'f and s 
and asymptotically at large s >> 4mQ it depends only on Let ns stress that the 
dependence of 11 on s is crncial at low initial energies only. 

The invariant inclusive spectrum can be also presented in the following equivalent 
form: 

1 da _ 1 da 
d^p TT dp'jdy TT drnl^dy 

Taking into account flTTl) we can rewrite Eq. flldl) in the form: 

n dmjdy ^ = 0, H) + (j)g{y = 0, H) ■ (1 - and/ 9 {{s/so)^)] ■ g{s/so)^ . (18) 

The first part of the inclusive spectrum (Soft QCD (quarks)) is related to the function 
4>q{y = 0) n), which is fitted by the following form [30] : 

0g(2/ = O,n) = Aqexp{-U/Cq) , (19) 

where Ag = 3.68 {GeV/c)~‘^,Cq = 0.147 GeV/(?. 

The function (j)g{y = 0,11) related to the second part (Soft QCD (gluons)) of the 
spectrum is fitted by the following form [3U] : 


(j)g{y = 0,U) = Ag^/m^texp{-Il/Gg) , 


( 20 ) 


where Ag = 1.7249 {GeVlc)-\Gg = 0.289 GeV/c\ 

Using flTSi) we can calculate the inclusive hadron spectrum as a function of the trans¬ 
verse mass. 

In Figs, m Ejthe inclusive spectra {1 /mit)da/dmudy of 7r“-mesons produced in pp col¬ 
lisions at the initial momenta Pin = 31 GeV/c and Pin = 158 GeV/c are presented versus 
their transverse mass rrit. Using only the first part of the spectrum ^^(y = 0,mit), which 
is due to the quark contribution, the conventional string model, let’s call it the SOFT 
QGD (quarks), one can describe the NA61 data [TT] rather satisfactorily at Pin =31 GeV/c 
and rrit < 1 GeV /G . This part of the inclusive spectrum corresponds to the dashed line 
in Fig. |2j The inclusion of the second part of spectrum due to the contribution of gluons 
(SOFT QGD (gluons)), the dotted line, allowed us to describe all the NA61 data up to 
rrit 1.5 Gev/c, see the solid line in Fig. |2] (Soft QGD(quarks-l-gluons)). Actually, at large 
^/s even at the NA61 energies If ~ mit/rriQ instead of 0161) . Generally the pion spec¬ 
trum p{s,mii = E{d^a/d^p) (ignoring the gluon part) can be presented in the following 
approximated form, which is valid for the NA61 energies and low transverse momenta 
Pt < 1 GeVja. 


p(s,miQ ~ (\)g{y = 0,n)^(g/go)^ = g{s/ sq)^A qexp{- ^ mo(l^^ 4mg/s) ^ 

= 9{s/sQ)^Aqexp{-^) , 

where 


T = Gqmo{l — 4mo/s) (22) 

is the inverse slope parameter, which is called sometime as the thermal freeze-out tempei- 
ature. One can see from Fq.([22|) that this thermal freeze-out tempeiaXme depends on the 
initial energy square s in the c.m.s. of collided protons. That is the direct consequence 
of the self-similarity approach, which operates the four-momentum velocity formalism. 
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P|n=31 GeV/c 



Figure 2: Results of the calculations of the inclusive cross section of hadron production in 
pp collisions as a function of the transverse mass at the initial momenta Pin = 31 GeV/c. 
They are compared to the NA61 experimental data HU. 


This s-dependence of T is signihcant at low initial energies and at s >> tuq the slope 
parameter T becomes independent of s. To describe rather well the NA61 data at larger 
values of pt, the inclusive pion spectrum should be presented by Eq. flT8|l . which has a more 
complicated form in comparison to the simple exponential one of (23). However, the main 
contribution to the inelastic total cross section comes from the hrst part of Eq. flT8|l . which 
has the form given by Eq.(21) 

We have calculated the inverse slope parameter T, as a function of the energy s/s given 
by Eq. fl22p and presented in Fig. 01 There is a good agreement with the experimental 
data III]. 

In Figs. [SEwe give the calculations of inclusive spectra of charged hadrons (mainly 
pions and kaons) produced in pp collision at ^/s = 900 GeV, 2.36 TeV, 7 TeV performed 
by using flTT]) and the perturbative QCD (PQCD) within the LO [SniEI] compared to the 
LHC data. These spectra are the sums of inclusive spectra of pions and kaons, therefore 
they are presented as a function of the transverse momentum pt instead of functions of 
the transverse mass mu because the masses of a pion and kaon are different. In addition 
to the part of spectrum, which corresponds to (IT^ . see the solid lines in these hgures, we 
also include the PQCD calculations, see the dotted lines. The PQCD calculation within 
the LO is divergent at low pt, therefore, the dotted lines go up, when pt decreases. The 
kinematical region about pt — 1.8-2.2 GeV/c^ can be treated as the matching region of 
the nonpertubative QCD (soft QCD) and the pertubative QCD (PQCD). One can see 
from Figs. [5]- [7] that it is possible to describe rather well these inclusive spectra in the 
wide region of pt at the LHC energies matching these two approaches. Figs. I71IH1IM 7I show 
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Pin=158 GeV/c 



Figure 3: Results of the calculations of the inclusive cross section of hadron production in 
pp collisions as a function of the transverse mass at the initial momenta Pin = 158 GeV/c. 
They are compared to the NA61 experimental data HU. 


that in order to describe the inclusive hadron spectra [I71|TS] in the mead-rapidity region 
at rriit < 2 GeV/c^, the form of spectra presented in a simple exponential form of ffTTD 
should be changed and can be presented in the form of flT^ . which includes the non trivial 
energy dependence. To describe rather well the LHC data on these inclusive pt-spectra 
at pt > 2-3 GeV/c, the PQCD calculation should be included, the contribution of which 
has a shape similar to the power law p^-distribution [ 25 ] . 


5 Conclusion 

The inclusive hadron spectrum in the space of four-velocities is presented within the 
self-similarity approach as a function of the similarity parameter H. The use of the 
self-similarity approach allows us to describe the ratio of the total yields of protons to 
anti-protons produced in A-A collisions as a function of the energy in the mid-rapidity 
region and a wide energy range from 10 GeV to a few TeV. To study the similar ratio of 
light nuclei to anti-nuclei, we need more detailed experimental data. 

We have shown that the energy dependence of the similarity parameter H included 
within this approach is very signihcant at low energies, namely at \fs < 6 GeV, and 
rather well reproduces the experimental data on the inverse slope or the thermal freeze- 
out temperature of the inclusive spectrum of hadrons produced in pp collisions, it increases 
and saturates when -s/s grows. This is very signihcant for a theoretical interpretation of 
future experimental data planned to get at FAIR, GBM (Darmstadt, Germany) and NIGA 









a> 

Q. 

O 


0 ) 

> 

c 


in 

a> 

V) 



pp -> n x - 

NA61 data at 0.0 < y < 


0.2 ^ 


2 


4 


6 8 10 12 14 16 18 20 



Figure 4: Results of calculations of the inverse slope parameter T on the energy depen¬ 
dence for the negative pion production in pp-interactions. The experimental points are 
taken from m. 


(Dubna, Russia) projects. That is an advantage of the self-similarity approach compared 
to other theoretical models. 

However, we have also shown that the s dependence of H is not enough to describe 
the inclusive spectra of hadrons produced in the mid-rapidity region, for example, in pp 
collisions in the wide region of initial energy, especially at the LHC energies. Therefore, we 
modify the self-similarity approach using the quark-gluon string model (QGSM) [261127] 
and [MU including the contribution of nonperturbative gluons, which are very signihcant 
to describe the experimental data on inclusive hadron spectra in the mid-rapidity region at 
the transverse momenta pt up to 2-3 GeV/c [HOlIM] . Moreover, the gluon density obtained 
in [31] , parameters of which were found from the best description of the LHC data, allowed 
us also to describe the HERA data on the proton structure functions [12] . To describe the 
data in the mid-rapidity region and values of pt up to 2-3 GeV/c, we modify the simple 
exponential form of the spectrum, as a function of H, and present it in two parts due to 
the contribution of quarks and gluons, each of them has different energy dependence. This 
energy dependence was obtained in [30| based on the Regge approach valid for soft hadron- 
nucleon processes. To extend the application of the suggested approach to analyze these 
inclusive pf-spectra at large hadron transverse momenta, we have to include the PQCD 
calculation, which results in the main contribution at pt > 2-3 GeV/c. 


References 


[1] E.Fermi, Phys. Rev. 92 (1953) 452. 

[2] 1. Ya. Pomeranchuk, Izv. Dokl. Akad. Nauk Ser.Fiz. 78 (1951) 889. 

[3] L.D. Landau, Izv. Akad. Nauk Ser. Fiz. 17 (1953) 51. 


9 







s^^^=900 GeV 



Figure 5: Results of the calculations of the inclusive cross section of charge hadrons 
produced in pp collisions at the LHC energies as a function of their transverse momentum 
Pt at y/s =0.9 TeV. The points are the LHC experimental data [T8] . 

[4] R. Hagedorn, Supplemento al Nuovo Cimento Volume HI, 147 (1965); 

R. Hagedorn, Nuovo Cimento 35 (1965) 395; 

R. Hagedorn, Thermodynamics of strong interactions,Nuovo Cimento 56 A (1968) 
1027 

[5] J. Cleymans, H. Oeschler, K. Redlich and S. Wheaton, Phys. Rev. C73 (2006) 034905; 
arXiv: 0511094 [hep-ph]. 

[6] A. Andronic, P. Braun-Munzinger, J. Stachel, Nucl. Phys. A 834 (2010) 237C; 
arXiv:0911.4931 [nucl-th]. 

[7] . A. Andronic, P. Braun-Munzinger, J. Stachel, Nucl. Phys. A 772 (2006) 167; 
arXiv:0511071 [nucl-th]. 

[8] F. Becattini, J. Manninen, and M. Gazdzicki, Phys. Rev. C 73 (2006) 044905; 
arXiv: 0511092 [hep-ph]. 

[9] S. V. Akkelin, P. Braun-Munzinger, Yu. M. Sinyukov, Nucl. Phys. A 710 (2002) 439; 
arXiv:0111050 [hep-ph]. 

[10] M.E.J. Newman, Contemporary Physics 46, 323 (2005). 

[11] A. Clause!, C.R. Shalizi, M.E.J. Newman, SIAM Review 51, 661 (2009). 


10 







s^^^=2.36 TeV 



Figure 6: Results of the calculations of the inclusive cross section of charge hadrons 
produced in pp collisions at the LHC energies as a function of their transverse momentum 
Pt at y/s =2.36 TeV. The points are the LHC experimental data |18] . 

[12] F. Bardou, J.-P. Bouchaud, A. Aspect, C. Cohen-Tannoudji, “Levy Statistics and 
Laser Cooling’, Cambridge University Press, (2002). 

[13] M. Mitzenmacher, Internet Mathematics, 1, 226 (2003). 

[14] B. 1. Abelev et al. (STAR Collaboration), Phys. Rev. C 75, 064901 (2007). 

[15] A. Adare et al. (PHENIX Collaboration), Phys. Rev. C 83, 052004, (2010); Phys. 
Rev. C 83, 064903 (2011). 

[16] K. Aamodt, et al. (ALICE Collaboration), Eur. Phys. J. C 71 1655 (2011); Phys. 
Lett. B693, 53 (2010) ; Phys. Rev. D 82, 052001 (2010). 

[17] G. Aad, et al. (ATLAS Collaboration), New J. Phys. 13, 053033 (2011). 

[18] V. Khachatryan, et al. (CMS Collaboration), Phys. Rev. Lett. 105, 022002 (2010). 

[19] 1. Bediaga, E.M.F. Curado, J.M. de Miranda, Physica A 286, 156 (2000). 

[20] G. Wilk and Z. Wlodarczyk, Eur. Phys. J. A40, 299 (2009); 48, 161 (2012). 

[21] R. Hagedorn, Riv. Nouvo Cimento 6, 1 (1984). 

[22] C. Michael and L. Vanryckeghem, J. Phys. G 3 L151 (1977). 

[23] C. Tsallis, J. Statist. Phys. 52, 479 (1988). 


11 







S^^^=7TeV 



Pt, [GeV/c] 


Figure 7: Results of the calculations of the inclusive cross section of charge hadrons 
produced in pp collisions at the LHC energies as a function of their transverse momentum 
Pt at ^/s =7 TeV. The points are the LHC experimental data pTlITS] . 


[24] C. Tsallis, R. S. Mendes, A. R. Plastino, Physica A 261, 534 (1998). 

[25] J. Cleymans, G.I. Lykasov, A.S. Parvan, A.S. Sorin, O.V. Teryaev, Phys.Lett. B723 
(2013) 351; arXiv: 1302.1970 [hep-ph]. 

[26] A. B. Kaidalov, Z.Phys. C12 (1982) 63. 

[27] A. B. Kaidalov, Surveys High Energy Phys. 13 (1999) 265. 

[28] B. Andersson, G. Gustafson, G. Ingelman, T. Sjstrand, Phys.Reports, 97 (1983) 31. 

[29] K.Werner, Phys.Reports, 232 (1999) 87. 

[30] V. A. Bednyakov, A. A. Grinyuk, G. 1. Lykasov, M. Pogosyan. Int.J.Mod.Phys., A27 
(2012) 1250042. 

[31] A. A. Grinyuk, G. 1. Lykasov, A. V. Lipatov, N. P. Zotov. Phys.Rev. D87 (2013) 
074017. 

[32] A. M. Baldin, L. A. Didenko. Fortsch.Phys. 38 (1990) 261-332. 

[33] A. M. Baldin, A. 1. Malakhov, and A. N. Sissakian. Physics of Particles and Nuclei, 
Vol.32. Suppl. 1, 2001, PP.S4-S30. 

[34] A. M. Baldin, A. 1. Malakhov. JINR Rapid Gommunications, 1 [87]-98 (1998) 5-12. 


12 







[35] A. M. Baldin, A. A. Baldin. Phys. Particles and Nuclei, 29 (3), (1998) 232. 

[36] A. Tawfik. Nuclear Physics A 859 (2011) 63-72. 

[37] http: //hepdata. cedar. ac. uk/view/p7907 . 

[38] R. Klingenberg et al. Nuclear Physics A 610 (1996) 306c-316c. 

[39] V. Abramovsky, V. N. Gribov and O. Konchelli, Sov. J.Nucl.Phys. 18 (1973) 308. 

[40] K. A. Ter-Martirosyan. Sov.J.Nucl.Phys., 44, 817 (1986). 

[41] A. A. Abgrall et al. Eur.Phys.J., C74 (2014) 2794. 

[42] A. V. Lipatov, G. I. Lykasov, N. P. Zotov. Phys.Rev. D89 (2014) 014001. 


13 



